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ABSTRACT 

We have investigated past measurements of the local supermassive black hole 
mass density, correcting for hitherto unknown dependencies on the Hubble constant, 
which, in some cases, had led to an underestimation of the mass density by factors 
of ~2. Correcting for this, we note that the majority of (but not all) past studies 
yield a local supermassive black hole mass density that is consistent with the range 
4.4-5.9 xlO 5 /( J ff o )M Mpc~ 3 (when using H = 70 km s" 1 Mpc -1 ). In addition, we 
address a number of ways in which these past estimates can be further developed. In 
particular, we tabulate realistic bulge-to-total flux ratios which can be used to estimate 
the luminosity of bulges and subsequently their central black hole masses. 
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1 INTRODUCTION 

Assuming that the dark mass concentrations at the centres 
of both elliptical galaxies and the bulges of disc galaxies 
are the sleeping engines that powered past quasar activity 
(e.g. Miller 2006; Brand et al. 2005, and references therein), 
then the local mass density of such quiescent supermassive 
black holes (SMBHs), Pbh,o, can be used to constrain mod- 
els of quasar formation and growth (e.g. Haehnelt & Kauff- 
mann 2001; Bromley, Sommerville & Fabian 2004; Yu & Lu 
2004; Hopkins, Richards & Hernquist 2006). After factoring 
in potential SMBH mass-energy losses due to gravitational 
radiation (Ciotti & van Albada 2001; Yu & Tremaine 2002; 
Menou & Haiman 2004) and the possibility of lost mass from 
either "three-body" SMBH slingshot ejection (e.g. Volonteri, 
Haardt & Madau 2003; Hoffman & Loeb 2007) or explusion 
via gravitational radiation recoil (e.g. Merritt et al. 2004; 
Libeskind et al. 2006) pbh.o helps constrain the amount of 
material to explain past quasar flux. 

Together with the quasar luminosity function (e.g. Hop- 
kins et al. 2005, and references therein) integrated over time, 
Pbh.o can also constrain the average efficiency at which mat- 
ter is converted to radiation as it falls onto a SMBH (e.g. 
Ciotti, Haiman & Ostriker 2001; Elvis et al. 2002; Ferrarese 
2002; Yu & Tremaine 2002; Fabian 2003; Marconi et al. 2004; 
Merloni 2004; Shankar et al. 2004; Yu & Lu 2004). This can 



in turn tell us about the rotation of SMBHs. For exam- 
ple, a non-rotating Schwarzschild black hole is expected to 
have an efficiency of 5.4 per cent while a maximally rotat- 
ing Kerr black hole may have an efficiency as great as 37 
per cent (Thorne 1974; Hasinger 2005). Radiative efficien- 
cies are typically reported to range around 10—15 per cent 
but values as high as ~ 30—37 per cent are also sometimes 
reported (e.g. Gallo et al. 2004; Crummy et al. 2006; Wang 
et al. 2006). 

For the above reasons it is of interest to accurately de- 
termine Pbh.o- In an attempt to help explain some of the 
differences between previously reported values (Table [TJ, we 
will discuss a number of corrections and adjustments that 
could be made to past estimates. We focus on how estimates 
of Pbh.o depend on the Hubble constant, and how past mea- 
surements which have not fully taken this into account are 
affected — sometimes changing by factors of 2 or more. We 
shall refer to these revised estimates of pbh,o (Section [2]) as 
our "/i-corrected values". In addition we raise a number of 
other points pertaining to the accurate estimation of the 
local SMBH mass density, mostly addressing the issue of 
recovering the host bulge luminosity before converting this 
into a SMBH mass. Section [3] provides a summary. 



2 ILLUSTRATIVE STUDIES 
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The nature of the hidden or over-looked dependencies on 
the Hubble constant are endemic to most past estimates of 
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Table 1. Local SMBH mass density estimates. The factor h 3 = [H /(70 km s _1 Mpc" 1 )] 3 
is appropriate for the Graham et al. (2007) study because the M^—n relation they used is 
independent of the Hubble constant. The majority of the densities from other papers have been 
transformed to Hq = 70 km s _1 Mpc -1 using h? rather than h 3 , as indicated in each paper. 
However, as shown in Table [3] this is not always appropriate. 



Study 


p b h,o (E/S0) 


Pbh,o (Sp) 


Pbh,o (total) 




h% 10 5 M Q Mpc" 3 


h^ 10 5 M Q Mpc" 3 


h^ 10 5 M e Mpc' 3 


Graham ct al. (2007) 


(3.46 ± 1.16)/i 7 o 


(0.95 ± 0.49)/i 7 o 


(4.41 ± 1.67)/i 7 


Wyithe (2006) 


(3.4+^)^ 




2.28 ± 0.44 


Fukugita & Peebles (2004^ 






Marconi et al. (2004) 


3.3 


1.3 


^• D -1.4 


Shankar et al. (2004^ 


q 1+0.9 

-0-8 


i 1+0.5 
± - 1 -0.5 


42 +l.l 


Shankar et al. (2004^ 


3 0+ 10 


1 2+ - 4 


42 +l.l 


McLurc & Dunlop (2004) 


2.8 ±0.4 




9 9 +3.9 


Wyithe & Loeb (2003) 






Aller & Richstone (2002)£ 
Yu k, Tremaine (2002)^ 


1.8 ±0.6 


0.6 ±0.5 


2.4 ±0.8 


2.0 ±0.2 


0.9 ±0.2 


2.9 ±0.4 


Merritt & Ferrarese (2001)£ 






4.6fef 1 


Salucci ct al. (1999) 


6.2 


2.0 


8.2 



a See their equation 75. 
^ Based on their Section 3.2. 
c Based on their Section 3.4. 
d Taken from their Table 2. 

e Based on their (M"bh — (j)-derived mass function. 
/ See also Ferrarese (2002). 

Pbh,o- It is therefore necessary to only look at a couple of 
representative case studies in detail, and provide the revised 
estimates from other studies in tabular form (Table [3} . In 
what follows, we have chosen two interesting and well writ- 
ten studies. 

2.1 Case study 1 

Our first example is the analysis by Aller & Richstone (2002, 
hereafter AR02), who used the L-a relation (e.g. Faber & 
Jackson 1976) to convert luminosities into velocity disper- 
sions and then applied the Mbh-ir relation (Ferrarese & Mer- 
ritt 2000; Gebhardt et al. 2000) to obtain a histogram of 
SMBH masses. We have identified two areas for improve- 
ment pertaining to the treatment of the Hubble constant. 

The first is in regard to their adjustment of the 
Tremaine et al. (2002) Ad^-a relation for what they referred 
to as a correction from h — 0.8 to h — 1, where h — f/o/(100 
km s -1 Mpc -1 ). Had all, or at least the majority of, the 
SMBH masses used to construct the Tremaine et al. rela- 
tion been obtained with distances that depended on an as- 
sumed Hubble constant of 80 km s -1 Mpc -1 , then it would 
be appropriate to multiply (decrease) the SMBH masses by 
a factor of 0.8 to make the M^-a relation consistent with a 
Hubble constant of 100 km s -1 Mpc -1 (AR02, their equa- 
tion 23). However, only five of the 31 galaxies used to con- 
struct the Tremaine et al. relation had distances, and thus 
SMBH masses, derived using a Hubble constant of 80 km 
s -1 Mpc -1 ; most galaxies had their distances obtained using 
surface brightness fluctuations (Tonry et al. 2001). Removal 
of the five galaxies with distances obtained using Hq — 80 
km s -1 Mpc -1 changes neither the slope nor intercept of 
the Tremaine et al. M^-a relation by more than 0.01. This 
relation is therefore effectively independent of the Hubble 



constant and need not be adjusted. Consequently, the total 
SMBH mass density in AR02 should be 25 per cent higher 
and scale with h 3 rather than h 2 . Similarly, the /i-correction 
in Yu & Tremaine (2002, after their equation 6) which was 
applied to Tremaine et al.'s (2002) M\^-a relation should 
also not have been applied. 

The SMBH masses that were computed by AR02 are 
dependent on their adopted Hubble constant for a second 
reason: their SMBH mass estimates were derived from abso- 
lute magnitudes which depend on h. This SMBH dependence 
on h can be seen in the x term which appears in their equa- 
tion 25, and which is defined in their equation 22. Removing 
the aforementioned factor 0.8/Zi from their equation 25, one 



has Mbh oc h 



-4.02x5/7.7 
-2.61 



oc h 



That is, their M* term 



varies with h~ . This has apparently gone over-looked in 
the literature to date. 

Now, AR02's equation 24 for the SMBH mass function, 
which is in units of Mbh (rather than log Mbh) and which 
depends on M„ , can be written as 



dN 



dMbh M t h 



h 3 4>* 



-2.61 



Mb 



-2.61 



exp 



f M bh y 



■(1) 



The expression for the SMBH mass density is thus 

Pbh 



Afbh-m— dMbh 



dMi 



bh 



-2.61 



P 
a + 2 



a±2 



M mi] n 



13 



t M min y 



/ M max 
}. 



(2) 



where y(a,x) is the incomplete gamma function (e.g. Press 
et al. 1992) defined by 
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7(a, x) = I e t t a dt. (3) 
Jo 

The SMBH mass density in AR02 therefore actually varies 
with ft 3 " 2 ' 61 = h - 39 (not h 2 as given in AR02) and it 
also varies with an additional complicated dependence on 
h which is tied up in the gamma functions above. Correct- 
ing AR02's SMBH mass density to H = 70 km s" 1 Mpc -1 , 
and integrating down to a minimum mass of 1O 6 M0, as they 
did, one obtains a value of Pbh = 5.9 x 10 5 M o Mpc -3 . This 
value is ~2.5 times larger than what AR02 find when h = 0.7 
(see their equation 28). 

This type of correction is again not unique to the anal- 
ysis in AR02, for example, McLure & Dunlop's (2004) esti- 
mate of pbh for E/SOs, increases from 2.8 xl0 5 M Q Mpc -3 
to 4.8xl0 5 M o Mpc" 3 (for h = 0.7, see Graham et al. 2007, 
their Section 4). In general, all SMBH mass functions which 
have been derived from /i-dependent galaxy luminosities will 
depend on h in a similar fashion, although the above factor 
of 2.61 may vary from paper to paper (see Table [3}. 



Table 2. Mean ± standard deviation of the (bulge minus galaxy) 
magnitude and bulge-to-galaxy flux ratios, B/T, derived from lit- 
erature data. (A more complete summary, using other literature 
data, and a variety of optical and near-infrared passbands will be 
presented in Graham & Worley 2007, in prep.) 

SO/SOa Sa,Sab,Sb >Sc 
R 1,/4 -bulgc plus Exponential-disc 



(mbuige - mtotH 0.61 ±0.32 1.37 ±0.68 3.22 ±0.99 

(B/T) a 0.59 ±0.16 0.34 ±0.20 0.07 ± 0.06 

(m bu i ge - mtotli ■■■ 1.67 ± 1.06 3.18 ±1.41 

(B/T) b ... 0.32 ±0.26 0.11 ±0.11 

Exponential-bulge plus Exponential-disc 

(m bu i go - mtot) 6 ... 2.70 ± 1.18 4.08 ± 1.03 

(B/T) b ... 0.12 ±0.10 0.04 ±0.03 

/j!/ ?l _bulge plus Exponential-disc 

(mbulgc - mtot)^, ... 2.36 ± 1.06 4.21 ± 1.06 

(B/T) c ... 0.17 ±0.09 0.03 ±0.03 

(B/T)£ 0.25 ±0.09 
( B / T )L 0.24 ±0.11 



2.1.1 Related issues 

Ignoring the above mentioned dependencies on the Hubble 
constant for the moment, it is expected that the SMBH 
masses in AR02 are too high at the low-mass end because of 
a) the way they converted disc galaxy magnitudes into bulge 
magnitudes, and b) the way they assigned a velocity disper- 
sion to these magnitudes. The average bulge-to-disc (B/D) 
luminosity ratios which AR02 assigned to their early- and 
late-type spiral galaxies, and also lenticular galaxies, came 
from the R 1 ' 4 -bulge plus exponential-disc decompositions in 
Simien & de Vaucouleurs (1986). Due to Simien & de Vau- 
couleurs use of the R 1 ^ model to describe bulges which are 
better matched with an R}^ n profile having n <~ 3, and 
often around 1 (e.g. Andredakis & Saunders 1994; de Jong 
1996; Balcells et al. 2003), too much flux has been assigned 
to the bulges of their disc galaxies. We have derived the 
mean bulge-to-total ratios from various studies and show 
the results in Table [2] On average, the bulge luminosities 
used by AR02 will be ~2 times too bright and thus their 
estimate of the SMBH mass in the bulges of disc galaxies 
will be high by a factor of ~2. Correcting for this would re- 
sult in a 12.5 per cent reduction to their total SMBH mass 
density, giving a value of 5.2 xlO 5 M Mpc -3 (# = 70 km 
s -1 Mpc -1 ). Using Allen et al.'s (2006) Sersic-bulge + ex- 
ponential disc decompositions of 10,095 galaxies, we plan to 
apply the Mbh-L relation from Graham (2007) to obtain a 
new measurement of p b h,o- 

Regarding the conversion of these overly-bright bulge 
magnitudes to velocity dispersions, the logarithmic slope of 
the L-a relation is known to be shallower at fainter luminosi- 
ties (e.g. Tonry 1981, Held et al. 1992), and for magnitudes 
below Mb ~ — 19.5 ±1 mag the slope is approximately 2 (de 
Rijcke et al. 2005; Matkovic & Guzman 2005), compared to 
a value of four for the more luminous spheroids (Faber & 
Jackson 1976; their figure 16). The use of a constant slope 
of 3 by AR02 would have therefore systematically under- 
estimated the velocity dispersion as one progresses to fainter 
magnitudes, and over-estimated the velocity dispersion in 
the larger spheroids. Without performing a full re-analysis 
of their data, the overall corrective term is unknown. In pass- 



a B-band data from Simien & de Vaucouleurs 1986. 
^ B-band data from de Jong 1996. 
c B-band data from Graham 2003. 

^ _fC-band data from Balcells, Graham & Peletier 2007. 
e X-band data from Laurikainen, Salo & Buta 2005. 

ing we note that the non-linear nature of the L-a relation 
complicates AR02's prediction of the parameter (3 shown in 
their equation 10. 

2.2 Case study 2 

Our second example is the analysis in Shankar et al. (2004, 
their section 3.1 & 3.2), who used the Mbh~L relation, in 
addition to the M^-a relation, to estimate pbh from var- 
ious luminosity functions. Their equation 1 for predicting 
SMBH masses from luminosities was also obtained under 
the (false) assumption that the SMBH masses which de- 
fine this relation are dependent on the Hubble constant. 
They modified the M\^-L relation from McLure & Dun- 
lop (2002; their equation 6) which had originally been (cor- 
rectly) constructed with no i/o-adjustment to the black hole 
masses. The equation in Shankar et al. therefore requires 
that log(70/50) be subtracted from the left hand sid^E Cor- 
recting this results in a 40 per cent increase to their (Mbh~ 
L)-estimated SMBH masses and thus a 40 per cent increase 
in their value of pbh- Their /i-corrected value for pbh, total is 
(5.9 ± 1.5) x 10 5 M o Mpc" 3 (/i = 0.7), in perfect agreement 
with AR02's /i-corrected value. Similarly, their equation 3 
should not contain the factor 80/i/o, and so their (M^-a)- 
estimated SMBH masses, and thus their (Mbh-cr)-derived 

1 Equation 1 from McLure & Dunlop (2004) should read 
log(Mbh/M Q ) = 1.251og(Lx/i0) - 5.53, when using M K & = 
3.28 mag and Ho = 70 km s -1 Mpc -1 . While relevant to the pbh 
adjustments made in Table [3] this is perhaps a moot point given 
that the galaxy distances and thus absolute magnitudes used to 
construct that equation are known independently of the Hubble 
constant, just as the SMBH masses are, see section 3.1 of Graham 
2007. 
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/5bh,o, needs to be increased by 14 per cent when using their 
adopted value ol h = 0.7. 

Shankar et al.'s equation 4 for the SMBH mass function, 
which has c6* in units of log Mbh (rather than Mbh) can be 
written as 



diV 



d log Mbh 



= h 3 <j>, ( 



M„h- 



exp 



/ Mbh \ ' 
VM./i" 2 - 5 / 



(4) 



where the exponent —2.5 comes from their equation 1 which 
was used to transform magnitudes Mr into SMBH masses 
using log A/bh oc —0.5(Mr + 5 log h). The expression for the 
SMBH mass density is thus 



Pbh 



h 3 <j} : 



-7 



M»h~ 



diV 



d log Mbh 
'a + 2 



d log Mbh 



/3(ln 10) 
a + 2 



( Mmax y 

Vm*/i- 2 - 5 J 



/3 



/ M min y 
Vm*/i- 2 - 5 / 



}■ 



(5) 



A similar parameterisation of the SMBH mass function 
could be made for the data in McLure & Dunlop (2004) 
and Marconi et al. (2004), except for the latter study the 
exponent would be —2.26 rather than —2.5 (see their equa- 
tion 10). This full dependency on the Hubble constant was 
not included in Tundo et al.'s (2007) reanalysis of these 
works. Their Mbh oc L 1,30 relation is also considerably 
steeper than the new expression Mbh oc L - 93 reported by 
Graham (2007) and it predicts notably larger SMBH masses 
for galaxies more luminous than Mr ~ — 21 mag. 

In Table [3] we provide updated values of pbh and, im- 
portantly, show their dependence on the Hubble constant. 
While some estimates of pbh.o appear not to have changed 
from Table [T] one should note that the quoted dependence 
on h may have changed, which is of course of importance if 
H 70 km s" 1 Mpc" 1 . 



2.2.1 Related issues 

Many studies have assumed the universal existence of R 1 ^ 
light-profiles when obtaining their total galaxy magnitudes, 
and have thus introduced a systematic bias into their 
luminosity-derived SMBH mass function. For instance, a 
light-profile shape dependent — and therefore luminosity 
dependent (e.g. Graham & Guzman 2003, their figure 10 and 
references therein) — magnitude correction (Graham et al. 
2005) is applicable to the SDSS Petrosian magnitudes which 
Shankar et al. used. Adding —0.2 mag to the Petrosian mag- 
nitudes, in an effort to recover the total galaxy magnitude, 
is only applicable if every galaxy has an R 1 ^ 4 light-profile. 
However, a range of profile shapes has long been known to 
exist (e.g. Davies et al. 1988; Caon et al. 1993) and is such 
that a smaller/greater correction for missed flux needs to be 
applied to the Petrosian magnitudes of galaxies less/more 
luminous than Mr ~ —21 mag (Kormendy & Djorgovski 
1989). Similarly, a light-profile shape (and outermost sam- 
pled radius) dependent magnitude correction (Graham & 
Driver 2005, their figure 10) is required for recovering total 
magnitudes from Kron magnitudes. Indeed, half a galaxy's 
flux may be missed using Kron magnitudes (Andreon 2002; 
Bernstein, Freedman, & Madore 2002; Benitez et al. 2004). 



3 SUMMARY 

Table [3] shows our "ft-corrected" pbh,o values. It should be 
noted that the ft-dependent corrections we have detailed ef- 
fect not only the value of Pbh.o but also the SMBH mass 
functions from which these values are typically derived. The 
related issues we have raised in sections 12. 1.1 1 and 12. 2. II that 
pertain to the luminosity of the host spheroid have not been 
included in Table [3] The use of ft- independent Sersic indices 
and velocity dispersions for constructing the SMBH mass 
function and mass density results in a purely h 3 dependence 
for pbh- This is because the SMBH masses that are involved 
are derived from relations which themselves do not depend 
on any assumed Hubble constant. 

We (tentatively) identify previously missed agreements 
on the value of pbh,o- For example, AR02's corrected valu^fl 
of (5.9 ± 2.0) x 10 5 M Q Mpc~ 3 (for h = 0.7) is now in 
good agreement with Merritt & Ferrarese's (2001) (Mbh~ 
cr)-derived measurement of 4.6 x 10 5 Ii 3 Mq Mpc -3 . In fact, 
a (near) consensus on the local SMBH mass density now 
exists. The M^-L based studies are seen to agree with 
each other and with recent studies which have used a 
mean SMBH-to-spheroid mass ratio convolved with the lo- 
cal spheroid mass density. The M^-n based study (Gra- 
ham et al. 2007) is also seen to agree with both of these 
types of analysis, with the optimal (total) SMBH mass den- 
sities ranging from 4.6-5.9 xlO 5 M Mpc~ 3 (/i = 0.7) for all 
three types of analysis. Furthermore, some of the ft-corrected 
Mbh-o" based studies (Aller & Richstone 2002; Marconi et 
al. 2004; Shankar et al. 2004) also provide consistent results 
with this range. The two exceptions are the noticeably lower 
values of (2.9 ± 0.4) x 10 5 /i^ M Q Mpc -3 (Yu & Tremaine 
2002) and (2.0 ± 0.4) x lQ 5 h 3 70 M & Mpc~ 3 (Wyithe 2006)0. 

Excluding galaxies without 'secure' SMBH mass de- 
terminations, Marconi et al. (2004) derived and used an 
Mbh-o" relation with a 0.17 dex higher zero-point (and 0.09 
steeper slope) than used by Yu & Tremaine (2002). This ac- 
counts for their different (Afbh-o")-derived values of pbh,o. It 
is also worth noting that if the local sample of ~30 galaxies 
with direct SMBH mass measurements have low luminosities 
with respect to the greater population at any given veloc- 
ity dispersion (Yu & Tremaine 2002; Bernardi et al. 2007; 
Tundo et al. 2007; Lauer et al. 2007), then the M hh -L rela- 
tion will over-predict pbh,o- However, as noted by Graham 
(2007, his Appendix), until accurate bulge/disc decomposi- 
tions are available for the greater population, and correc- 
tions for dust attenuation in the bulges of disc galaxies are 
addressed (Driver et al. 2007), this remains uncertain. 
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2 Computing AR02's disc galaxy's bulge luminosities with more 
realistic B/T flux ratios leads to a 12.5 per cent reduction in 
Pbh,o. giving a value of (5.2±1.8) x 10 5 M Q Mpc" 3 (for h = 0.7). 

3 Correcting Wyithe's (2006) estimate of pbh f° r the intrinsic 
scatter in the M^-a relation (see equation 12 in Yu &; Tremaine 
2002) would only increase it by a factor of 1.1-1.2. 
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Table 3. Modification of Tabic [T] Here, the local SMBH mass density estimates have been fully corrected for their 
dependence on h, and transformed to Ho = 70 km s — 1 Mpc - 1 . While some estimates of pbh,o appear not to have 
changed from Table [l] one should note that the quoted dependence on h may have changed. The term f(h) is used 
to denote that a more complicated dependence on h exists and needs to be taken into account if one is to transform 
these values to a different Hubble constant (see, e.g. equation [2}- 



Study 


Method 


Pbh.o (E/S0) 


Pbh.o (Sp) 


Pbh.o (total) 








lO 5 A/ Mpc" 3 


10 5 A/ Q Mpc" 3 


1O 5 M Mpc" 3 


Graham ct al. (2007) 


A*bh- 


-n 


(3.46 ± 1.16)h 3 


(0.95 ± 0.49)h 3 


(4.41 ± 1.67)/i 3 


Wyithe (2006) 


M bh - 


(7 


(1.98 ± 0.38)ftf 


Fukugita & Peebles (2004) 


Psphcroid 


(3A+_\ i - ! )h 70 


(1.7i^)ft 70 


(5.1+ 3 -*)h 70 


Marconi et al. (2004) 


M bh -(L, a) 


3.3h%™f(h) 


1.3h% r4 f(h) 




Shankar et al. (2004) 


M hh - 


L 


(4.3±{-l)h° 5 f(h) 


(l.5±° 7 7 )h° Q Sf(h) 


(5.9±[l)h%Sf(h) 


Shankar et al. (2004) 


M bh - 


a 


(3-4_ 7 )/i 3 




(4-81^)^0 


McLure & Dunlop (2004) 


M hh - 


L 


(4.8 ± 0.7)/i° 6 5 / (/i) 




(2.1^)^ 


Wyithe & Locb (2003) 


M bh - 


(7 






Aller & Richstone (2002) 


M bh - 


(7 


(4.5 ±1.5)h° 39 f(h) 


(1.4±1.3)h° 39 /(h) 


(5.9 ± 2.0)h%™f(h) 


Yu & Tremaine (2002) 


M bb - 


a 


(2.0 ± 0.2)/i 3 


(0.9 ± 0.2)/i 3 Q 


(2.9 ± 0.4)/i 3 


Merritt & Ferrarese (2001) 


Psphcroid 






4.6/170 


Salucci et al. (1999) 


Psphcroid 


G.2h% n 


2.0h% a 


8.2/i| n 



Appendix of a different paper. We also thank A.Marconi, 
D. Merritt, F. Shankar and Q.Yu for their comments. 
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